Heat cost of parametric generation of microwave squeezed states 
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*5 ! Abstract 

In parametric systems, squeezed states of radiation can be generated via extra work done by 
external sources. This eventually increases the entropy of the system despite the fact that squeezing 
is reversible. We investigate the entropy increase due to squeezing and show that it is quadratic 
in the squeezing rate and may become important in the repeated operation of tunable oscillators 
(quantum buses) used to connect qubits in various proposed schemes for quantum computing. 
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Parametric devices using the nonlinear inductance and natural protection from deco- 
lerence of Josephson junctions are considered as prospective elements of quantum circuits 



Jj. In particular, these have been proposed and implemented as quantum buses, providing 
tunable coupling between solid state qubits (see, e.g., (2-6]). On the other hand, para- 
metric Josephson devices can produce squeezed states of microwave radiation UQ. Such 



experiments were performed in the 1980s 12j, [l3| and more recently 



0. 



In this paper we investigate the generation of squeezed states during the operation of a 
quantum bus and estimate their additional contribution to the entropy of the system and 
associated heat production. We will see that this contribution is quadratic in the squeezing 
rate and should be taken into account if the system operates near its limit of efficiency 



15 



16|). 



determined by the Landauer erasure principle (see, e.g. 

A quantum bus can be considered as a harmonic oscillator with tunable frequency oj(t); 
tuning it in and out of resonance with qubits allows to manipulate their states (e.g., by 
performing two-qubit gates). The changes of u(t) should be fast compared to the charac- 
teristic frequencies of qubits and the bus, i.e., u/u ^> u; the decoherence rate in the circuit, 
however, must be much smaller than u. But these are the very conditions which lead to the 
generation of squeezed states. Indeed, a coherent state of a harmonic oscillator is squeezed 



by a sudden change of the oscillator frequency [17] . Recall that a coherent state is a state 



with equal and minimal uncertainties: (AP 2 ) = (AQ 2 ) = 1/2, where P 



[a — a 



)/iy/2 



and Q — (a + aJ)/ a/2 are dimensionless momentum and position expressed through bosonic 
creation and annihilation operators. 



It is convenient to describe oscillator states by their Wigner functions (see, e.g., jl8|): 

1 



W(a,a*) 



2vr 2 



dX dX* e (- Aa * +A * a ) tr 



e (\at-\*a) p 



(1) 



which have the advantage of reducing to classical distribution functions in the classical limit. 
Here p is the density matrix of the system. For a coherent state, the Wigner function is 
Gaussian, 



W ao (a,a* 



— exp[— 2\a — ao\ 



7T 



The amplitude and phase of the complex parameter «o (which in the Schrodinger represen- 
tation behaves as a (t) = a (0) exp[— iuit]) describe the classical limit of the oscillator state. 
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A squeezed state will have instead the Wigner function 

2 f 2 2 2 1 

W(a, a*) = — exp < —2s [(x — xq) cos 9 + (y — yo) sin 9] \{y — y ) cos 9 — (x — xq) sin 9} > 

7T I s J 

(2) 

Here x = Re a, y = Im a, s is the squeezing parameter, and 9 determines the direction of 
the squeezing axis. Similarly, the thermal state 

2/ TV 



l + 2n 



r exp[-2|a|7(l +2fi)], 



characterized by the average photon number n = {exp [uj/T] — 1} 1 , can be squeezed to fl9^ 
2/tt 



W(a,a* 



l + 2n 



exp 



l + 2n 



s (x cos 9 + y sin 9) 2 H — (y cos 9 — x sin 9) 2 



(3) 



After being introduced in Ref. [17J, squeezing the oscillator states 



oscillator frequenc y w as considered in a number of papers 



was shown 
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jy a sudden change of 



2CH24j. In particular, it 



24( that, in the absence of decoherence, repeated abrupt small changes of 



the oscillator frequency can produce arbitrarily large squeezing. Here we will mainly follow 



the approach of 



Consider an arbitrary time dependence of the system frequency u(t). Let us also denote 
the creation/annihilation operators belonging to the Fock state of an oscillator with the 
frequency u(t = 0) by a , aj,, and those corresponding to uj(t), by a u , a^. The Hamiltonian 
keeps its standard form, 

H(t) = hu(t) (4a w + i) , (4) 

and the commutation relations between a w , a\ hold, if the old and new operators are related 
via a Bogoliubov transformation 

[u(t) + w(0)] a - [u(t) - w(0)] aj 



2y/u(t)u(0) 



(5) 



which can be explicitly written as [1 



a = V*{t)a u V{t); aj = V\t)alV{t) 



(6) 



where 

V(t) = exp 



In 



exp - 



In 



(7) 
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In order to work in the initial Fock space at t = 0, we apply the transformation V 1 (t) 
V'(t), which gives 

H(t) -> H{t) = V(t)H(t)V\t) - iM{t)^-V\t) = 



hu{t) ( ala + - ) + ifi—j-r cLq — (al)' 2 ■ (8) 



2/ £j(t) 

Hereafter we correct an error made in 



t\2 



11], where the transformation (JSJ) was effectively 
applied twice in the same direction, which quantitatively (but not qualitatively) affected the 
results. For the Wigner function, with a, a* always refer ring to the coherent states in the 



same Fock space (at t = 0), we find the master equation ( ll|, mutatis mutandis) 



^W(a,a*,t) = 2w(f)Bn W(a,a*,t) + ^^Re (a-^) W(a,a*,t), (9) 

or 

9 Tir , \ ,J d d\ Tjr/ , ld\nu(t) { d d\ Tjr/ , 

«" , (*.».«) = «W W'^ + i-ar 1 ^- "Ty) w{x ' y ' t] - (10) 

We did not include the diffusive terms, which describe decoherence (including relaxation), 
thus assuming that the system maintains coherence over many cycles of frequency- switching 
and of its free evolution with frequency u(t). Equation ffTO]) is a first-order linear equation 
and can be solved by the method of characteristics. The characteristic equations are 

dx Co dy Co , . 

— = —x -uy; — = ux- —y. 11 
dt Zoo dt 2uo 

In the limit of fast frequency change, when the cu-terms dominate, these equations lead to 
squeezing of the quantum state: 

d In x din JZo d\ny d In Juj , , 
~ i . z. ~ v (12) 

dt ~ dt ' dt ~ dt ' V ; 

so that W(x, y, At) ~ W(y/sx, -^y,0), with squeezing parameter s = u(At)/uj(0), and At 
the duration of the fast frequency change [cf. Eqs. (|2J [3}]. In the quasistatic limit, Eq. ( ITT|) 
simply describes the rotation of the Wigner function as a whole: slow changes of parameters 
do not produce any squeezing, as expected. 

As we mentioned above, in the absence of decoherence, periodic changes of the oscillator 
frequency can produce an arbitrarily high degree of squeezing, even if the difference betweeen 
the two limiting values, oo and u>i, is arbitrarily small. It is only necessary that at least 



one of the transitions (cuq — > uj\ or ujq — > oj\) is fast on the scale of u, and these transitions 
are tuned to the phase of the oscillator. This may also take place during the operation of 
a quantum bus, which would normally switch fast on- and off-resonance with the qubits 
coupled to it. As we shall see, there are entropy costs associated with the production of 
squeezed states. Therefore this side-effect should be avoided in the operation of quantum 
buses. 

Let us return to the characteristics equation ffTTj) . Without loss of generality, we can 
consider two regimes of oscillator frequency switches uq -h- uj\\ the "ratchet" regime (when 



one transition, e.g., oj —tux, is fast, and the opposite one is slow [11|, |24(), and the "seesaw" 
regime (when both transitions are fast 21] ). In either case, if two consecutive fast switches 
occur at the moments t n _i and t n , the point on a characteristic will evolve according to 
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(13) 



tn-1+0 



Here Q n ,n~\ is the phase angle accumulated during the period of slow evolution, and s n is 
the squeezing achieved at the nth step. In the "ratchet" regime: s n = s = oj\/ojq (or vice 
versa), while in the "seesaw" case: s 2n = l/s 2n +i = s - Obviously, detA n = 1. After two 
consecutive switches 



[S - 1) COS n , n _l COS 9 n -l >n -2 + cos 6 l „ in _2 

i - 1J cos0 n , n _i sin n _i,„_ 2 - sin0 n)n _ 2 
and 



s - 1) sin^ nin _i sin + COS 6> n , jn _2 

i - l) sin6> i cos6> n _i n _ 2 - sm0 nin _ 2 



(A n A„_!) r = 
s - 1) cos# nin _i sin + Sin 6 n ,n-2 

- 1) COS 9 n .„_i COS n _i n _2 + COS n , n _2 



(14) 



(A n A n _O s 

[s - 1) sin n , n _i cos n -i,n-2 + sin 0„, n _2 



1) sin0„ n _i sin0 n _i re _ 2 + cos6> _ 2 



.(15) 



Here n , n _ 2 
to 



@n,n-i + @n-i,n-2- In the case of periodic switchings, n +i )Tl = 9, this reduces 



(A„A 



n-l, 



's-1) cos 2 6 + cos 29 



^±sin29 



l/s+l 



sin 29 



i- l") cos 2 9 + cos 29 



(16) 



and 



(A n A n _ 



- 1) sin 2 9 + cos 29 



Vs+l 



sin 29 



^±1 sin 29 
(7-1) sin 2 9 + cos 29 



(17) 



The resonance conditions, O r = nq and s = n(q+l/2), when either (1T6|) or ( TT7|) diagonalize, 
ensure that after 2N switchings the state will be squeezed exponentially to an arbitrary 
degree, 





± 



/ s N x 



s~ N yo 



(18) 



no matter how small (s — 1). For finite values of (s — 1), an exact resonance is not necessary. 



For example, in the ratchet case a runaway squeezing will happen if 24 1 

1 + I sin 61 ,„ _ 

| cos 6 1 

Deviations from periodic switching in the operation of a quantum bus, fluctuations of circuit 
parameters, and eventually relaxation and decoherence in the system will limit the actual 
degree of squeezing. Nevertheless, the very operation of a quantum bus will tend to produce 
squeezed states. 

Squeezing is in itself a reversible process. Nevertheless, in the presence of decoherence 
it will lead to an increase of the system's entropy, which can be associated with internal 



friction |25] . To be specific, consider the so-called "energy entropy" 

Se[p\ = -J2p^ ln P« ^ s Ip\ = _tr \p ln p] ' 
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26] 



(20) 



where p is the density matrix, {p n }^Lo are its diagonal terms in the energy basis, and S is 
the standard (fine-grained) von Neumann entropy. The two entropies coincide if and only if 
p commutes with the Hamiltonian. Obviously, unlike S, Se can increase in a closed system. 
Decoherence processes eventually wipe out the off-diagonal components of the density matrix 
in the energy representation, thus increasing the von Neumann entropy to the level of Se- 
This makes Se a convenient tool for the study of non-equilibrium states, especially since 
unlike the von Neumann entropy, Se can be conveniently expressed through the Wigner 
function [27(, since 



Pn " 

2(-l) n 



= 2(-l) n J da da* e" 2|Q|2 L n {4\a\ 2 ) W{a,a*) 
J J dxdye' 2 ^^ L n (A(x 2 + y 2 ))W(x,y). 



(21) 



Here L n (x) is the Laguerre polynomial. We can associate with the additional entropy, 
SS(s) = Se(s) — Se(0), of a squeezed state a specific amount of heat, 



SQ(s) = T 6S E (s), 
6 



(22) 



which will be eventually released into the system. 

As an example, consider a squeezed thermal state ([3]), for which an exact solution can be 



actually found 27J in terms of the hypergeometric function 



28 



= «£<*<-!>■ (t^)" +1 " 8 .F, (i,„ + 1 - l; -=£^£>) , (23) 

where k = tanh(w/2T). A more convenient approximation, valid for small squeezing, k(s — 
1) <C 1, is given by 

"-^(»(i^ + T^))-^(tC + ^). <-) 

Here = (1— exp[— oj/T\) exp[— nu/T] is the equilibrium population of the nth energy level 
of the oscillator, e = \s — 1| characterizes squeezing, Iq(z) is the modified Bessel function, 
and the last expression in Eq. (I24p is valid for oj/T <C 1. The normalization constants 
Ci,2 = 1 + 0(e 2 ). The occupation of states with numbers n < Nq will initially decrease, 
and of those with numbers n > Nq increase with squeezing. Expanding the Bessel function, 
I ( z ) = 1 + z 2 /4 + . . ., we see that d(p n - p c ^)/de = 0, if 

2T 1 

„ = N^---. (25) 

This is in good agreement with the results based on the exact formula in Eq. ( 123]) (Fig. [I]). 
The corresponding change in entropy is given by 

(co/T) 2 



5S E (e) = -J2$Pnlnp: 



eq _ 

16(1 -exp[-u;/r]) 2 



E 2 [ , UJ S , CJ UJ ] f UJ 



-|(l-- )ln --3 + 5-j [-«!). (26) 
With the same accuracy, the equilibrium entropy of an unsqueezed thermal state becomes 

S = - ln[2 sinh(w/2T)] + (u/2T) coth(uj/2T) « - ln(w/T) + I. 
Therefore in the leading term in T/s ^> 1 the squeezing contribution to the entropy is 

5S E (e) ~ jS, (27) 

and the amount of additional heat released into the system due to squeezing will be 

e 2 (s - I) 2 
5Q ~ T £ -S = { -^-^TS (28) 
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FIG. 1: (Color online.) Occupation numbers of the energy levels of a harmonic oscillator in a 
squeezed thermal state [Eq. ([23]) ] as a function of squeezing e = |s — 1|. (a) w/T = 0.2, estimated 
boundary number N = 2T/uj - 1/2 = 9.5; levels (top to bottom) 8, 9, 10, 11, 12. (b) oj/T = 0.4, 
iVo = 4.5; levels (top to bottom) 3, 4, 5, 6, 7. 



S 



per cycle. Detuning between the bus and qubits can be typically 10-30% {gj, making the 
prefactors in Eqs. ( 127|28|) as high as 0.25, i.e., the contribution of squeezing to the system 
entropy becomes comparable to the equilibrium entropy of the device. Of course, this 
contribution becomes important only when the system operates near the Landauer erasure 
limit, with Ql = Tin 2 15|, |l6|. If a need ever arises to limit its effect, the quadratic 



dependence of 5Q on the squeezing parameter indicates that even a slight decrease of the 
difference between the operating frequencies of the bus will drastically reduce it. One 
also should avoid a periodic operation of the bus, to exclude the resonant increase of the 
squeezing rate. 



CONCLUSION 



We have shown that the operation of a quantum bus will lead to squeezing of its 
quantum state, and calculated the corresponding additional heat, which will be injected 
into the system. This contribution is quadratic in the squeezing rate (proportional to the 
ratio of the working frequencies of the bus) and can become important only in systems 
operating either near the Landauer erasure limit or under the conditions when there is a 
resonant increase of squeezing by repeated bus operation. 
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